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The electronic and magnetic properties of the cubic double perovskite Ba2NaOsOe were examined 
by performing first-principles density functional theory calculations and analyzing spin-orbit coupled 
states of an Os 7+ (d 1 ) ion at an octahedral crystal field. The insulating behavior of Ba2NaOs06 was 
shown to originate from a cooperative effect of electron correlation and spin-orbit coupling. This 
cooperative effect is responsible not only for the absence of orbital ordering in Ba2NaOsOe but also 
for a small magnetic moment and a weak magnetic anisotropy in Ba2NaOsOs. 

PACS numbers: 71.70.Ej, 71.27. +a, 71.20.-b, 75.30.Gw 



Oxides of orbitally degenerate 3d transition metal ions 
at octahedral sites exhibit rich electronic and magnetic 
properties arising from the interplay between spin, or- 
bital and charge degrees of freedom. In an octahedral 
crystal field, the d-orbitals of a transition metal ion are 
split into the t2 g and e g levels. A 3d perovskite with un- 
evenly filled e g levels, e.g., LaMnC>3 with high-spin Mn 3+ 
(3d 4 ) ions, has a strong tendency for orbital ordering []]]. 
In contrast, a 3d perovskite with unevenly filled ti g levels, 
e.g., YTi03 with Ti 3+ (3d 1 ) ions, has a reduced tendency 
for orbital ordering and hence provides opportunities to 
observe an intricate interplay between the spin and or- 
bital dynamics @, Q • It is an important issue to under- 
stand the mechanisms that select the ground state out of 
numerous possible states arising from this competition 

dsaii. 

Most studies probing this question have focused on 
3d oxides, and much less is known about whether re- 
lated 4d and 5d oxides can exhibit similar behavior. 
Due to a large spatial extension of 5d orbitals, effects 
of electron-correlation are weaker in 5d oxides than in 
3d oxides. However, effects of spin-orbit coupling (SOC) 
are stronger in 5d oxides than in 3d oxides. Thus, 5d 
oxides should exhibit a different balance between spin, 
orbital and charge degrees of freedom. In this context, 
it is of interest to examine the electrical and magnetic 
properties of the cubic double perovskite Ba2NaOs06 
@, 0, Si- In this 5d oxide the NaC>6 octahedra share 
corners with the OsC>6 octahedra of orbitally degenerate 
Os 7+ (5d*) ions, and the nearest-neighbor OsC>6 octahe- 
dra run along the [110] direction while the next-nearest- 
neighbor OsC>6 octahedra run along the [100] direction 
(Figure [T|). Ba2NaOsOg presents several puzzling prop- 
erties. It is an insulator despite the fact that the struc- 
ture remains cubic down to 5 K without any distortion of 
the OsOe octahedra from their regular octahedral shape 
[8] . It is unclear what mechanism lifts the orbital degen- 
eracy of the Os 7+ (d 1 ) ions to make Ba2NaOs06 insu- 
lating. The magnetic susceptibility of Ba2NaOs06 be- 
tween 75 and 200 K follows a Curie- Weiss law with a 
negative Weiss temperature (i.e. 9 as —10 K), which 



TABLE I: Comparison of the zero-field magnetic moments 
l-iexp (per FU) with the calculated moments jj, ca i c (per 
FU) for the ferromagnetic state of Ba2NaOs06 using the 
GGA+SOC+U method with U eff = 0.2 Ryd a - b . 
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0.19 


[110] 


0.51 


-0.34 


0.98 


0.64 


0.22 


[100] 


0.52 


-0.30 


0.97 


0.67 


0.18 



a /is(Os) and /it(Os) are the spin and orbital magnetic 
moments calculated for the Os atom, respectively, and 
Us is the calculated total spin magnetic moment per FU. 
b The moments are in units of /is. 
c Ref. [8] 

shows that the dominant spin exchange interaction be- 
tween Os 7+ ions is antiferromagnetic (AFM) 0, [|| . How- 
ever, Ba2NaOs06 undergoes a ferromagnetic (FM) order- 
ing below Tq = 6.8 K with a very low magnetic moment, 




FIG. 1: (Color online) Structure of cubic double perovskite 
Ba 2 NaOs0 6 . 

To gain insight into the puzzling electronic and mag- 
netic properties of Ba2NaOs06, we carried out a first- 
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principles density functional theory (DFT) electronic 
structure study. In this Letter, we show that the in- 
sulating property of Ba 2 NaOs0 6 originates from a co- 
operative effect of electron correlation and SOC. This 
cooperative effect is responsible not only for the lack of 
structural distortion in Ba 2 NaOs06 but also for the low 
magnetic moment and the weak magnetic anisotropy Q 
of Ba 2 NaOs0 6 in the FM state. 

Our first principles DFT electronic structure calcu- 
lations were performed by using the full-potential aug- 
mented plane waves plus local orbital method as imple- 
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FIG. 2: (Color online) Band structures calculated for 
Ba 2 NaOs0 6 using different methods, (a) GGA, (b) GGA+U 
with U eff = 0.2 Ryd, (c) GGA+SOC, (d) GGA+SOC+U 
with U e ff = 0.2 Ryd. In (a) and (b) the solid and dashed 
lines refer to the up- and down spin bands, respectively. In 
(d) the valence band top is taken as the zero-energy point. 



mented in the WIEN2k code [9(. The non-overlapping 
muffin-tin sphere radii of 2.50, 2.14, 1.86, and 1.65 au are 
used for the Ba, Na, Os and O atoms, respectively. The 
expansion in spherical harmonics of the radial wave func- 
tions were taken up to Z = 10. The value of R™/^K max 
was set to 7.0. The total Brillouin zone was sampled with 
125 k-points. For the exchange-correlation energy func- 
tional, the generalized gradient approximation (GGA) by 
Perdew, Burke and Ernzerhof [101 ] was employed. The 
SOC was included on the basis of the second-variational 
method using scalar relativistic wave functions [11| . The 
structural parameters of Ba 2 NaOsOg were taken from 
the experimental values [Tj . 

The spin-polarized GGA band structure calculated for 
the FM state of Ba 2 NaOs06, presented in Figure [^a), 
has the Fermi level crossing both the up- and down- 
spin i 2g bands. This does not agree with the fact that 
Ba 2 NaOs06 is an insulator Q. At T the three i 2g bands 
are degenerate. The overall width of the i 2g bands is nar- 
row (approximately 1 eV) because the nearest-neighbor 
O- • O distance between adjacent OsOg octahedra is long 
(3.216 A) compared with the van der Waals distance of 
3.04 A. The exchange splitting of Ba 2 NaOsOg is about 
0.34 eV, which is considerably smaller than typical values 
found for 3d magnetic oxides. Because the bandwidth is 
small compared with the on-site repulsion U (approxi- 
mately 3.3 eV), it was suggested [8j that electron cor- 
relation is important in Ba 2 NaOs06- The failure of the 
traditional DFT in describing strongly correlated systems 
is currently remedied by the DFT plus on-site repulsion 
U method [13, E3]- Thus, we employed the GGA+U 
method to see if the insulating property of Ba 2 NaOs06 
can be explained. To avoid a double counting in the non- 
spherical part of potential, we used U e ff = U — J and 
omit the multipolar terms proportional to J in the added 
GGA+U potential. The band structure calculated for 
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FIG. 3: (Color online) DOS of Ba 2 NaOs0 6 obtained from 
the GGA+SOC+U calculation. The inset shows the PDOS 
plots calculated for the Os and O atom contributions to the 
t2 g bands. 
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the FM state of Ba 2 NaOs0 6 using the GGA+U method 
with U e ff — 0.2 Ryd is shown in Figure [2jb) , which re- 
veals a larger exchange splitting (about 0.7 eV). How- 
ever, there is no band splitting at T, and the dispersion 
characteristics of the ti g bands are almost the same as 
those of the GGA calculation. Our GGA+U calculations 
with larger U e ff values (up to 0.5 Ryd) did not change 
the general picture described above. Since Os is a heavy 
element, SOC is expected to play an important role in 
Ba2NaOs06- To see the effect of SOC on the electronic 
structure of Ba2NaOs06, we performed GGA+SOC cal- 
culations. Figure [2^c) shows the band structure calcu- 
lated for the FM state using the GGA+SOC method 
with the spin quantization taken along the [111] direc- 
tion. The up- and down-spin ti g bands are both split 
into three non-degenerate bands at T. However, these 
bands overlap with each other leading to a metallic state 
for Ba 2 NaOs0 6 . 

As described above, the GGA, GGA+U, and 
GGA+SOC methods all fail to reproduce the insulating 
state for Ba 2 NaOsOg. Nevertheless, we note that elec- 
tron correlation enhances the exchange splitting, while 
SOC splits the ti g bands. This suggests that a com- 
bined effect of electron correlation and SOC might in- 
duce a band gap in both the up- and down-spin t2 g bands. 
Therefore, we carried out GGA+SOC+U calculations for 
Ba2NaOsOg. The band structure calculated for the FM 
state with the spin quantization along the [111] direction 
is presented in Figure[2£d), which shows that Ba 2 NaOs06 
has an insulating gap with the lowest-lying down-spin i 2g 
band lying above the Fermi level. The associated density 
of states (DOS) calculated for Ba 2 NaOs06 is shown in 
Figure H The partial DOS (PDOS) plots for the Os 5d 
and O 2p states reveal that the Os 5d and the O 2p states 
contribute almost equally in the ti g bands. Though not 
shown, the PDOS plots calculated for the Os 5d xz , 5d yz 
and 5d xy orbitals show that these orbitals contribute 
equally to the t 2 g bands. The Os-0 bonding bands of 
the e g -symmetry occur well below the Fermi level (i.e., 
in the energy region between — 7 and —5 eV), which in- 
dicates the presence of a strong covalent bonding in the 
Os-0 bonds. 

So far, our calculations were performed for the FM 
state of Ba 2 NaOsOg. To find if the FM state is the mag- 
netic ground state, we considered an A-type AFM state 
in which the spins have the FM ordering within each 
sheet of Os 7+ ions parallel to the (001) plane, but have 
the AFM ordering between adjacent sheets parallel to the 
(001) plane. Our GGA+SOC+U calculation shows that 
the A-type AFM state is considerably less stable than the 
FM state (by 267 meV/FU), which is consistent with the 
observation that Ba 2 NaOs06 undergoes an FM ordering 
below 6.8 K j|. In the remainder of this work we will 
consider only the electronic structures calculated for the 
FM state of Ba 2 NaOs0 6 . 

To examine the effect of the spin quantization direc- 



tion on electronic structure, we carried out GGA+SOC 
and GGA+SOC+U calculations for the FM state of 
Ba 2 NaOs06 with the spin quantization taken along the 
[111], [110] and [100] directions. The stability depen- 
dence of the FM state on the spin quantization direc- 
tion is negligible in the GGA+SOC calculations, but is 
not negligible in the GGA+SOC+U calculations. In the 
latter calculations with U e ff — 0.2 Ryd, the FM states 
with the [110] and [100] quantizations are less stable than 
that with the [111] quantization by 1.9 and 16 meV per 
FU, respectively. The spin, orbital and total moments 
calculated for the FM state of Ba 2 NaOsOg using the 
GGA+SOC+U method with U eff = 0.2 Ryd are sum- 
marized in Table I. For each of the [111], [110] and [100] 
quantizations, the total moment per FU is calculated to 
be ~ 0.65/is, which is smaller than the spin-only value of 
1 hb due to the fact that the orbital moment of ~ 0.33/is 
is opposite to the spin moment in direction. Neverthe- 
less, the calculated values are still large compared with 
the zero-field moments (~ 0.2/ie/FU) determined from 
the magnetization study of Ba 2 NaOs06 0, 0] • 

The band gap opening at the Fermi level, the moment 
reduction and the slight magnetocrystalline anisotropy 
in Ba 2 NaOs06 can be accounted for by considering the 
effect of SOC on the i 2g orbitals of an Os 7+ ion. In the 
second variational approach for spin-orbit coupling 111 ], 
the scalar-relativistic part of the Hamiltonian is diagonal- 
ized on a basis adopted for each of the spin projections 
separately, and then the full Hamiltonian matrix is con- 
structed on the basis of the eigenfunctions obtained in 
the first step. The spin-orbit part of the Hamiltonian in 
the Os spheres is then given by 



H so — AL • S, 



(1) 



where the SOC constant A > for the Os 7+ (d 1 ) ion 
with less than half-filled t 2g levels. With 9 and 4> as the 
azimuthal and polar angles of the magnetization in the 
rectangular crystal coordinate system, the L • S term is 
rewritten as 

L-S = S Z (L Z cos d + \L + e-^ sin0 + sin 6) 



+ \S+{-L z s\xi0 - L+e 
+ \S-{-L z sin6> - L + e 



•i<j> 



sin 



+ L-e l f cos^f) 



cos 2 | + L_e^ sin 2 f). 



(2) 

Since the up- and down-spin i 2g bands are separated due 
to the exchange splitting, one can neglect interactions 
between the up- and down-spin states under the SOC to 
a first order approximation. This allows one to consider 
only the up-spin i 2ff bands using the degenerate perturba- 
tion theory, which requires the construction of the matrix 
elements (i\H so \j) = d xyi d yz ,d xz ) [Hj]. In this case, 
only the operators of the first line of Eq. [2] give rise to 
nonzero matrix elements. In evaluating these matrix el- 
ements, it is convenient to rewrite the angular parts of 
the d xy , d yz and d xz orbitals in terms of the spherical 
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harmonics as 







Y 2 - 2 ) 
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Using these functions, the matrix representation 
(i\H so \j) is found as 



iHX/2 



( 
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- sin 9 sin < 
sin 9 cos d> 



sm tt sin ( 


— cos# 




(4) 



Upon diagonalizing this matrix, we obtain the eigen- 
values of the three spin-orbit coupled states, namely, 
Ei = —hX/2, E 2 = 0, E 3 = h\/2. The associated eigen- 
functions ^1,^2,^3 are given by 



— [sin 9d xy + (i sin <f> — cos 9 cos 4>)d yz 



2 '"'""^s 
— (i cos<p + cos 6 sin (f>)d xz ] 

\&2 = ^[sm 9d xy — (i sin 4> + cos 9 cos (f))d yz 



(5) 



+(i cos — cos 9 sin ^)d az ] 

cos #eL + sin 9 cos 4>d yz + sin 9 sin < 



For these three states ^>2, and "3/3, the orbital mo- 
ments L along the spin quantization direction are — 1/J-b, 
and I^b, respectively, according to their eigenvalues 
and the SOC operator H so = AL • S given that S = 
1/2. Since the spin-orbit coupled state Vf'i is occupied, 
we obtain a negative orbital moment. The moment of 
— 0.35/Lts from the GGA+SOC+U calculation is consid- 
erably smaller in magnitude than — 1/xg. This is not 
surprising because the above analysis neglected the fact 
that the t 2g orbitals of an Os06 octahedron are not pure 
5d xz , 5d yz and 5d xy orbitals of the Os atom, but are 
given by their linear combinations with the 2p orbitals of 
the surrounding O atoms. Given that the Os 5d and the 
O 2p states contribute almost equally in the t2 g bands 
(Figure [3]), an orbital moment of approximately — 0.5/ib 
should be expected from the occupation of the up-spin 
t 2g band associated with v&i. 

The above analysis indicates that SOC splits the t 2g 
bands into three subbands regardless of the direction of 
the magnetization. However, SOC does not affect the 
exchange splitting. It is the on-site repulsion that en- 
hances the exchange splitting and increases the energy 
separation between filled and empty bands within each 
spin channel. That is, a cooperative effect of electron 
correlation and SOC is essential in opening a band gap 
at the Fermi level for Ba2NaOs06- Eq. [5] shows that 
the contributions of the d xy , d yz and d xz orbitals to the 
spin-orbit coupled state ^1 depend on the direction of 
the magnetization (i.e., 9 = 90° and </> = 0° for [100]; 
9 = 90° and cj> = 45° for [110]; 9 = arccos(V3/3) and 
4> = 45° for [111]). This is responsible for the weak mag- 
netic anisotropy observed for Ba2NaOs06 @] • The three 



d-orbitals contribute equally to the state *f?i for the [111] 
spin quantization, but unequally for the [110] and [100] 
spin quantizations. The inter-octahedron hopping inte- 
gral is nonzero along the directions of nearest-neighbor 
Os06 octahedra (e.g., [Ill] and [110]) but is practically 
zero along the directions of next-nearest-neighbor OsOe 
octahedra (e.g., [100]). For the [100] spin quantization, 
the ^1 level does not provide hopping along the [011] di- 
rection because it has no d yz orbital contribution. Con- 
sequently, the spin-orbit coupled band associated with 
should be higher in energy for the [100] quantization 
than for the [111] and [110] quantizations. This in part 
explains why the [100] quantization leads to a higher elec- 
tronic energy than do the [111] and [110] quantizations. 

In summary, the insulating behavior of Ba2NaOsOg is 
caused by a novel cooperative effect of electron correla- 
tion and SOC, which opens a band gap at the Fermi level 
and hence removes a driving force for orbital ordering of 
the orbitally degenerate Os 7+ (d 1 ) ions. The small mag- 
netic moment of Ba2NaOs06 arises from the fact that 
the occupied spin-orbit-coupled up-spin t 2g band gives 
rise to an orbital moment that is in opposite direction to 
the spin moment. 
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